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It has been recently proposed that the presence of a temporal electromagnetic field on cosmological 
scales could explain the phase of accelerated expansion that the universe is currently undergoing. 
The field contributes as a cosmological constant and therefore, the homogeneous cosmology produced 
by such a model is exactly the same as that of ACDM . However, unlike a cosmological constant term, 
electromagnetic fields can acquire perturbations which in principle could affect CMB anisotropies 
and structure formation. In this work, we study the evolution of inhomogeneous scalar perturbations 
in this model. We show that provided the initial electromagnetic fluctuations generated during 
inflation are small, the model is perfectly compatible with both CMB and large scale structure 
observations at the same level of accuracy as ACDM. 

PACS numbers: 98.80.-k, 95.36.+X, 14.70.Bh,04.62.+v 

I. INTRODUCTION 

Since the accelerated expansion of the universe was discovered more than a decade ago [l| , a big effort has been 
made in trying to identify the culprit responsible for it. The simplest explanation is provided by the introduction of a 
cosmological constant term in the Einstein equations. In fact, such a solution is in excellent agreement with most of the 
cosmological observations to date and only requires the introduction of one extra parameter to describe the energy 
density associated to the cosmological constant. Hence, as a phenomenological model to describe the observable 
universe, it is very successful. However, the value of the cosmological constant inferred from the observations is 
extremely tiny as compared to the other scale present in the gravitational action and set by the Newton's constant, 
namely the Planck scale. This poses a problem from the theoretical point of view because a theory containing two 
dimensional constants which differ by so many orders of magnitude does not seem to be very natural. This problem is 
usually referred to as the naturalness problem or the cosmological constant problem, and it is related to the coincidence 
problem, that is, why dark energy and matter have similar density parameters today when their respective energy 
densities have evolved very differently throughout the expansion history of the universe. These two problems are the 
main reason why so many models to explain the current acceleration of the universe have been proposed. Among 
those models, one can generally find models in which the accelerated expansion is caused by either a new field (dark 
energy) or a modification of Einstein gravity. The popular quintessence models 2] with a cosmological scalar field 
playing the role of dark energy, and their generalizations with non-canonical kinetic terms, like the k-essence models 
|3| or with a coupling to curvature invariants (see @ for a review on dark energy models) they all belong to the 
first group. In the second group we have the very popular f(R) theories @ or the models with extra-dimensions 0]. 
Although all these models were motivated by the cosmological constant problem, the actual situation is that most of 
them fail in solving it because they rather shift it to some other parameter of the theory. Besides, they are usually 
plagued by either classical or quantum instabilities or present conflicts with local gravity tests Q. 

Although the scalar fields have become the most popular candidates, vector fields have also been shown to be 
compelling candidates for dark energy 0, [Toil and even solve the coincidence problem [TT| . Generalisations to higher 
forms have also been considered recently Such models have been overlooked since they generically induce an 

anisotropy. This picture has changed in the last few years with the recent detections of some unexpected features 
in the CMB temperature anisotropies: A hemispherical as ym metry has been reported [l3l |: the angular correlation 
spectrum seems to be lac king power at the largest scales [Tj]; the alignment of the quadrupole and octupole (the 



so called Axis of Evil [15|, [16[ ) could also be seem as an extra-ordinary and unlikely result of statistically isotropic 
perturbations, even without taking into account that these multipoles happen also to be aligned to some extent with 
the dipole and with the equinox. The significance of the anomalies has been debated extensively in the literature (see 
e.g. 0, GJ, GJ, HJ EH, HI [H HI]) with some reported effects more significant than others [IE El, E3- 

One model that, not only succeeds in addressing the cosmological constant problem, but also is free from both 
classical and quantum instabilities and has the same set of PPN parameters as General Relativity is the one recently 
proposed in [28j |. In that model, the role of the cosmological constant is played by a temporal component of the 
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electromagnetic field 1 or, in other words, by the existence of an absolute cosmological electric potential [30(. That 
way, in addition of solving the aforementioned problems, we can establish the true nature of dark energy without 
resorting to unknown physics. In such a model, a slightly modified theory of the electromagnetic interaction is 
proposed and whose motivation comes from the quantization of the electromagnetic field in the covariant formalism. 
In that formalism, one modifies the electromagnetic action by introducing a gauge fixing term in the original action 
in exchange for either imposing the Lorenz condition d v A" — at the classical level or reducing the corresponding 
Hilbert space with the weak Lorenz condition d u A„ \4>) = at the quantum level. However, when moving to an 
expanding universe, the Lorenz condition cannot be consistently imposed at all times and, therefore, this approach 
becomes inappropriate. For this reason, in the model proposed in [28|, the electromagnetic interaction is quantized 
without imposing the Lorenz condition. The price to pay is that we need to introduce a new electromagnetic scalar 
mode in the theory so that the electromagnetic interaction would actually contain three physical degrees of freedom 
(see [U HH for other interesting consequences in theories with extra electromagnetic polarizations) . This new scalar 
mode is however decoupled from the conserved currents and can only be excited by gravitational fields. Then, it 
could be generated during an inflationary phase in the early universe and the interesting feature of the model is that, 
once it becomes super-Hubble, it gives rise to a new contribution to the action which is like that of a cosmological 
constant. Surprisingly enough is the fact that the observed value of the cosmological constant is naturally achieved if 
inflation occurred at the electroweak scale. 

Even though this model provides the same expansion history for the universe as a pure cosmological constant, its 
different nature makes possible to have a different evolution for the inhomogeneities so that it becomes necessary to 
perform an analysis of the perturbations to confront it to CMB and LSS data. In addition, as the observational data 
seem to be consistent with a cosmological term, it is then interesting that, by studying fluctuations, one in principle is 
able to infer the physical origin of the cosmological term. Previously the possibility of anisotropic fluctuations about 
a constant density term have been contemplated in the frameworks of viscous fluids and noncommutative geometry 
floL [29j . Here the properties of such a cosmological term can be consistently derived from a fully covariant theory. 
The presence of fluctuations can be understood since the constancy of the cosmological term is just a property of a 
solution for a component of a vector field. Uncovering the evolution of the fluctuations is the aim of the present work, 
which is organized in the following way. In Section I we give a short description of the model. In Section II we write 
down the background equations and solve them. Besides, we show how the background evolution gives rise to the 
same expansion history as a cosmological constant. Section III is devoted to the perturbations for the electromagnetic 
dark energy model. In particular, we obtain the corresponding equations and perturbed expressions in the Newtonian 
gauge and solve them in some interesting cases. The analogous expressions for the synchronous gauge are given in 
an appendix. Finally, in section IV we use a modified version of the publicly available numerical code CAMB (3f| 
to compute the evolution of the linear perturbations in the model. The results are summarized and discussed in the 
final Section IVll 



The proposed action describing the electromagnetic interaction including the gauge fixing term and current coupling 
is given by: 



with F^u — d^A u — d v A^. Here, the gauge fixing term is considered as a fundamental piece of the action on equal 
footing to the Maxwell term, i.e. electromagnetism is considered as a gauge non-invariant theory which contains three 
physical degrees of freedom, namely the two usual transverse photons plus an extra scalar mode (despite the gauge 
non-invariance, the theory is perfectly consistent as shown in [281]). In fact, the gauge fixing term can be seen as the 
kinetic term of this new mode and so the parameter £ can be fixed to 1/3 to have canonically normalized fields. While 
the general U(l) transformations are no longer an exact symmetry of the theory, it still preserves a residual gauge 
symmetry given by A^ — > A^ + d^9 provided O0 = 0. 
The modified Maxwell equations deduced from {1} read: 



II. DESCRIPTION OF THE MODEL 




(1) 



(2) 



1 We refer here to the vector as the electromagnetic field, instead of calling A M the vector potential and F^ v the photon. 
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If we now take the 4-divergence of these equations we obtain: 

n(y v A") = o 



(3) 



where we have used the fact that the electromagnetic current is covariantly conserved. The latter equation shows 
that V V A V behaves as a massless scalar field which is decoupled from the conserved electromagnetic currents, but it is 
non-conformally coupled to gravity so that it can be excited by gravitational fields. Moreover, due to the well-known 
fact that a massless scalar field gets frozen on super-Hubble scales for a FLRW universe, the four divergence of the 
electromagnetic field will be constant on scales larger than the Hubble radius, giving rise to a cosmological constant- 
like term in the action. On the other hand, for small scales we have that V „A V ~ a -1 so that its contribution to 
the Maxwell equations becomes negligible. Therefore, as long as the four divergence of the electromagnetic field is 
super-Hubble, it can play the role of a cosmological constant and, once it enters the horizon, it is rapidly suppressed 
and we recover the usual Maxwell equations with the Lorenz condition. 

It is interesting to note that the divergence of the electromagnetic field can be seen, at the equations of motion level, 
as a conserved current acting as a source of the usual Maxwell field. To see this, we can write — ^V^V^A") 
which, according to (|3|), satisfies the conservation equation V M J^. A = and we can express ((2|) as: 



On 



(4) 



with 



and V^Jy 



0. Physically, it means that, while the new scalar mode can only be excited by 
means of gravitational fields, once it is produced it will generally be considered as a source of electromagnetic fields. 
The energy-momentum tensor corresponding to the modified action {!]) is the sum of two pieces: 



T — T 

J-av — 



M 



fJ,IS 



(5) 



with 



T M — -n F aF a/3 — F F a 
the standard energy-momentum tensor corresponding to Maxwell theory and 

{V a A a f + 2A a V a (VpA )] - 4A^W u) (W a A a ] 



9p.v 



(0) 



(7) 



the energy- momentum tensor corresponding to the gauge fixing term. Notice that, for a residual gauge mode of the 
form A^ = d^9 with D0 = 0, the energy-momentum tensor vanishes identically. 



III. BACKGROUND EQUATIONS 



We shall assume a homogeneous and isotropic background described by the FLRW metric which, in conformal time, 
is given by: 

ds 2 = a(r) 2 [dr 2 - 8 l3 dx l dx3] . (8) 

In order to preserve the homogeneity we need the electromagnetic field to depend only on time whereas the large 
scale isotropy requires the absence of spatial components so that the background electromagnetic field is merely given 
by a temporal component which, in conformal time coordinates, we shall call Aq(t) 2 . Moreover, this homogeneous 
temporal component evolves according to the zero component of the equations @, that yield: 

A'a + 2 (H' - 2H 2 ) A Q = (9) 

with ' = 4- and 7i = — the Hubble expansion rate in conformal time. Notice that, since F^ v = for the assumed 
background configuration, the latter equation simply expresses that d(V „A V ) / dr = which, in turn implies that the 
four divergence of the electromagnetic field is constant irrespectively of the evolution of the scale factor. 



2 We shall denote by calligraphic letters the components of the vector field in conformal time coordinates 
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If we use now the fact that = together with V V A V = const in the expression for the energy momentum tensor, 
we obtain that: 



V = |(V^) 2 <^ (10) 

which is exactly the same as that of a cosmological constant of value A = | (V V A 1 ') 2 . Notice that in order to have 
a positive cosmological constant we need £ > 0. The background evolution for this model thus mimics that of the 
standard ACDM model whose cosmological constant value is determined by the four divergence of the electromagnetic 
field. The energy density of the electromagnetic field in this configuration is: 

PAo = ^ (A' + 2HA f (11) 

and the equation of state is simply pa = —PA - F° r the subsequent calculations, it will be convenient to introduce a 
redefinition of the vector field as Aq = a 2 Ao so that expressions (|9]) and (|ll|l look much simpler: 

(a- 4 A' Q y = (12) 
P*o = -^(A' ) 2 (13) 

and we clearly see that the energy density is constant. 

Finally, the fact that the background vector field does not contain spatial components is justified because on 
super-Hubble scales the spatial components grow more slowly than the temporal one. Let us see this in more detail. 
For a configuration with spatial components of the electromagnetic field in addition to the temporal one, i.e., with 
A^ = (Ao(t), A(t)) and in a FLRW universe, the modified Maxwell equations read: 

A% + 2(H' - 2H 2 )A = (14) 
A" = (15) 

Thus, if we assume a power law expansion with TL = p/t, the solution of these equations is: 

A = C iT 1+p + C 2T p (16) 
A = C x t + C 2 (17) 

and we see that the spatial components decay with respect to the temporal one in both the radiation (p = 1) and 
matter (p = 2) dominated eras. In addition, the energy density associated to the spatial part decays as a -4 whereas 
that associated to the temporal component is constant, as we have already seen. Hence, it makes sense to assume the 
aforementioned isotropic background with vanishing spatial components given above. 



IV. SCALAR PERTURBATION EQUATIONS 



In this section we shall derive the equations for the perturbations and solve them analytically for some simple cases. 
We only consider the case of scalar perturbations because we have checked that the vector perturbations evolve in the 
same way as in standard Maxwell theory and the tensor perturbations remain unaffected by the presence of the gauge 
fixing term as well. One might expect this because the gauge fixing term only affects the new scalar mode which, 
being scalar, can contribute to the scalar perturbations, but not to the pure transverse vector perturbations (i.e., the 
usual photons). We shall do the calculations of this section in the conformal Newtonian gauge, although we give the 
corresponding expressions for the synchronous gauge in the appendix. 

In principle, the effect of the high-electric conductivity of the universe should be taken into account including the 
corresponding electromagnetic current in the r.h.s of Maxwell equations. The reason why we can neglect such effects 
is the following: the electromagnetic current should satisfy the conservation equation: 

V M J" = (18) 

and also the condition of electric neutrality, i.e. 



(19) 
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with Ufj, the four-velocity of the comoving observers. Let us expand also the current up to first order as: 

J, = J° + 6J„ (20) 

but J® = for the homogeneous and isotropic electrically neutral background (this is the reason why we did not 
consider the current term in the background equations). If we now assume that the universe remains neutral at first 
order in the perturbations we have that 8 Jo = and, finally, from current conservation, we get: V • 8 J = 0, i.e. the 
perturbed current is transverse. In other words, when computing scalar perturbations, we can ignore the effect of 
electric conductivity by assuming that the electric charge vanishes at first order as well. 
In the considered gauge, the perturbed line element is given by: 

ds 2 = a(r) 2 [(1 + 20)dr 2 - (1 - 2(/))8 ij dx i dx j ] (21) 

In the absence of any anisotropic stress sources we have <fi = i>- In order to simplify the expressions we will use latter, 
we define the scalar perturbation of the vector field by introducing a convenient power of the scale factor as we did 
for the background field: 

6A„ = a- 2 (SA ,V6A) (22) 

As usual, we shall introduce the Fourier components of the perturbations and solve the corresponding equations for 
them. Then, from ^ we can obtain the following equations for the Fourier modes of the vector field perturbations: 



8A% k - m8A' ok - 3k 2 8A ok = ~2k 2 (28A' k - 5H8A k ) + M ~ W k + 30',' - A + 3(^ + cf>' k )A' ] 

(23) 



SA' k - mSA' k +(m 2 - 2H' - h 2 J SA k = 2 (^8A' ok - HSA ok 



1 \ 2 

3^ + ^41 A + -i> k A' 



In these equations we can see that the two perturbations of the vector field are coupled to each other and that 
the metric perturbations act as a source of them. Then, even in the case that the initial perturbations vanish the 
gravitational potentials will generate perturbations of the vector field that, eventually, may also be source of the 
metric perturbations. 

On the other hand, the corresponding perturbed energy-momentum tensor components are given by: 

ST °o = A {[- 2 V^o + SA' ok - (30' fe + ^' k )A ]A' + k 2 [(3SA' k - QH8A k - 3SA> k )A + 8A k A ]} 

STi j = A { l-^kA'o + 5A[ )k - (30; + ^' k )A ]A' + k 2 [(-3SA' k + 6H5A k + 36A ok )A + SA k A' Q ] } 6^ 

5T% = -^A Q [-(30i + 1>' k )Aa - 2^ k A' + SA' ak + k 2 SA k ] (24) 

It is interesting to note that this model has vanishing shear, i.e., 8T l - = —5p k <5*. This is due to there being only 
one physical scalar mode present. Moreover, from expressions (|24p one can find the following relation between the 
perturbed energy density and pressure of the field: 

2Ank 2 

5(pk + Pfc) = j- (SA ok - SA' k + 2H6A k ) (25) 

This relation is important because for a gauge mode satisfying a^ 2 SA^ = d^x (remember here the definition of the 
perturbations given in [j2"2"])) one has that 5Ao = 8 A' + 2H8A and, as a consequence, 8(p k + p k ) = 0. 

In order to obtain some analytical results, in the following we shall consider that the metric perturbations are 
generated by some dominating fluid and study the evolution of the electromagnetic perturbations in such a scenario. In 
other words, we shall assume that the perturbations of the electromagnetic field will not affect the metric perturbations 
evolutions. This assumption is justified as long as the electromagnetic energy density is clearly subdominant as it 
happens in most of the universe evolution when the energy density associated to the electromagnetic field is many 
orders of magnitude below that of the dominant component. However, such a condition will eventually breakdown at 
low redshift when dark energy becomes dominant and the results obtained here lack validity, being necessary to resort 
to a numerical treatment. In the early universe when radiation represents the dominant contribution to the energy 
density of the universe and neglecting neutrinos shear (which implies that ip = <f) the metric perturbation evolve as 

a ' 

Cife [lot cos(wt) - sin(wr)] + C 2k [ojt sin(wr) + cos(wt)] 
<Pk = 3 (2b) 
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Vector field perturbation in the radiation era Energy-momentum perturbation in the radiation era 




0.001 0.01 0.1 1 10 100 1000 0.001 0.01 0.1 1 10 100 1000 

kr kr 



FIG. 1: In these figures we show the evolution of the vector field perturbations in the Newtonian gauge for radiation (upper 
panel) and matter (lower panel) dominated universes. In the left panels we show the evolution of the vector field perturbations 
SAo/Ao (solid-blue) and kSA/Ao (dashed-red). The right panels show the perturbed energy density 5Tq (solid), pressure — \STl 
(dotted) and momentum — -^ki&T® (dashed). 



with lo = k/i/3. On the other hand, in the matter dominated universe, the gravitational potential becomes constant 
in time, i.e., 4>k = const. Then, we can solve the equations (|23|) in the presence of the gravitational perturbations 
produced by a radiation or matter fluid and obtain the evolution of the vector field perturbations in those epochs. 
The results are shown in Fig. [TJ We can see that, on super-Hubble scales, the perturbation 8Aq evolves in the same 
way as the background vector field so that SAq/Aq — const as one would expect. This also implies that the perturbed 
energy density is constant on large-scales as we can see in the figure. Moreover, this feature does not depend on 
the dominating fluid, i.e., it happens for both the radiation and matter eras. On small scales, the perturbed energy 
density scales with constant amplitude when the universe is dominated by radiation whereas the amplitude decays as 
1 /t in the matter era. Notice that this behavior on small scales is a common feature for all the perturbed components 
of the energy- momentum tensor, i.e., the energy, the pressure and the scalar component of the momentum. 

We shall obtain the explicit evolution for the electromagnetic perturbations during the matter-dominated era, when 
most of the cosmologically relevant scales reenter the horizon. In that epoch, the gravitational potential is constant 
as we said above so that = ipk — </>o and the Hubble parameter satisfies H' = — \J{ 2 . With these conditions, we 
can obtain the following expression for 8 A in terms of SAq and 4>q\ 



3 

oAk = ~TT? 



Mo* - Y nSA 'ok + (l k2 - 8 ^ 2 ) s <k - lnk 2 sA 0k - !4,<£o 



(27) 



This relation allows to find the following fourth-order differential equation for SAok'- 



sa 1 ^ - msAg + 2 ( k 2 - ^-n 2 ) SA'^ k - m 2 (k + 7 -n\ sA' ok + k 2 lk 2 + ^n 2 ) SA ok = 4k 2 nA' 4> 1 2s 1 



This equation together with the relation (|27p determine the evolution of the electromagnetic perturbations in a matter 
dominated universe. It is convenient to remind here the assumptions under which such equations remain valid. In 
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order to obtain those equations we have assumed that the metric perturbations act as an external source for the 
electromagnetic perturbations and that this external source is uniquely determined by the matter fluid. This means 
that the contribution of the electromagnetic field to the perturbed Einstein equations are negligible with respect to 
that of the matter component, which is a good approximation as long as the electromagnetic field energy density is well 
below the matter energy density. However, this condition eventually breakdowns because the electromagnetic field 
becomes dominant and it contributes in a non-negligible way to the Einstein equations so that the metric perturbations 
becomes affected by the electromagnetic field perturbations and the full system of coupled equations must be solved. 

In order to solve eq. (|28p we shall take advantage of the residual gauge symmetry of the theory, namely, — > 
A^ + 0^9 with HQ — 0. From this symmetry, we know that 5Aq — a 2 9' will be solution of (|28|) so that we can factorize 
it as .F[<?(<L4o)] = ik 2 HA' <j) , where J 7 is a second order differential operator and Q is the operator determining the 



evolution of a 2 9' which can be deduced from the equation satisfied for 9 and turns out to be Q 
Therefore, the equation (|28p can be expressed as: 



(It 2 



2H 



(It 



— -6H— 

dr 2 dr 



21 



dr 2 



d 



SAok = 4k 2 HA' <f> 



k 2 . 



(29) 



Thus, the solution for SAq will be determined by the equation Q(6Aq) = S with S the solution of the equation 
!F(S) = 4k 2 7i.A' (f)o- Since 7i = — in the matter era, we can obtain the explicit form of the kernel of T (solutions of 
the homogeneous equation) and is given by: 



S hom (r) = r 6 (C S ie- lkT + C S 2e lkT ) 



whereas a particular one can be obtained by: 

S 2 



S, 



part 



(30) 



(31) 



where Si and S2 are the two independent solutions given above and det Ws = S[S2 — S' 2 Si = 2ikr 12 CsiCs2 is the 
determinant of the Wronskian. Following the same procedure, we obtain the kernel of Q 



5A' ( 1 °" 1 {t) = d (k 2 r 2 - 3ikr - 3) e 
whereas the particular solution is 



-ikr 



C 2 {k 2 r 2 +3ikr - 3)e 



/7,-7 



8Al art = SA , Cl J 



SA 



C 2 



det W c 



S(r)dT - SA0.C2 



SAq, 



c 



det W c 



(32) 



(33) 



with 8Ao,Ci (* = 1) 2) denoting the piece of the homogeneous solution (f3"2")) proportional to C< (i = 1, 2) and det Wc — 
2ik 5 T C1C2. The homogeneous solution corresponds to the pure gauge degree of freedom and it does not contribute to 
the perturbed energy density, whereas the particular solution gives rise to the oscillating behavior with an amplitude 
decaying as r _1 for small scales shown in Fig. [1] 

We would like to remark that the procedure followed in order to solve the equations for a matter dominated universe 
is completely general and can be applied in any other situation where the metric perturbations is originated by some 
other component. 



V. EVOLUTION OF THE PERTURBATIONS 



In this section we shall present the results that we obtain from a modified version of the publicly available CAMB 
code [35{ to compute the CMB power spectrum when electromagnetic perturbations are taken into account. Since the 
background of the electromagnetic model is the same as that of ACDM, we do not need to modify the background 
equations of the code, although we do have to add the evolution equation for Aq . We use the first order equation (TT3|) 
so that we can relate the background electromagnetic field directly to the present value of its density parameter. The 
initial condition for this equation is unimportant for the background evolution, although it becomes relevant for the 
perturbations and it is set by assuming a power law behavior for Aq. For the perturbation equations, we add the two 
evolution equations for the electromagnetic perturbations given by J50|) and modify the corresponding terms involving 
dark energy perturbations. With these modifications, the code is ready to compute the evolution of the perturbations 
in the cosmology corresponding to the electromagnetic model and, thus, obtain the CMB power spectrum as well as 
the matter power spectrum. 
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Before proceeding to show the obtained results, we shall discuss what the initial conditions for the electromagnetic 
field perturbations should be. A natural origin for the presence of the new mode of the electromagnetic field on 
cosmological scales has been proved to be the quantum fluctuations of such a mode during an inflationary epoch 28] . 
In such a scenario, only the new scalar mode can be excited because of the conformal invariance of the usual transverse 
modes. The "homogeneous part" of this scalar mode, defined as the sum of all the modes which remain super-Hubble 
today, gives rise to the effective cosmological constant whereas the modes which have already reentered into the 
horizon constitute the origin of the electromagnetic perturbations discussed in the present work. In other words, 
we can split the primordial quantum fluctuations of the scalar mode generated during inflation in a homogeneous 
part comprising all the modes with k < ko (being kg the scale that is entering into the horizon today, i.e., the 
present Hubble radius) and an inhomogeneous part formed by those modes with k > k . Notice that such a split 
can be performed because the primordial power spectrum for the scalar mode is red-tilted, as shown in [28j | . so that 
the homogeneous part is large as compared to the inhomogeneous one and this enable us to treat the latter as a 
perturbation. For the mentioned scalar mode, one can see that the longitudinal component decays with respect to 
the temporal component on super-Hubble scales for a de-Sitter inflationary epoch (in an analogous manner to that 
shown for the background evolution) so that, at the end of inflation, the amplitude of the longitudinal component 
would be expected to be much smaller than the temporal one. Moreover, we have already shown that the longitudinal 
components also decay with respect to the temporal one in the radiation dominated epoch so that the super-Hubble 
modes would be expected to be strongly suppressed at the time when the initial conditions are given, which justifies 
to set the initial condition for SAk to zero. On the other hand, the power spectrum of the quantum fluctuations 
generated during a de Sitter inflation for the temporal component happens to be scale-invariant 3 so that we can set 
the initial condition for SAok as SAok (Tini) = A k~ z / 2 with A a constant depending on the details of the inflationary 
epoch such as the initial amplitude of the power spectrum after inflation or the duration of inflation. For our purposes 
in this work, this constant A will play the role of a free parameter to be constrained by comparing the CMB power 
spectrum produced by the model to the WMAP data. In the same way as for the background vector field, we shall 
give the initial condition for the derivative of SAok by assuming a power-law behavior. This is justified because we 
know that this is the type of evolution for the perturbations on super-Hubble scales, where the initial conditions are 
given. 

In Fig. [2] and Fig. [3] we show the results obtained from the modified version of the CAMB code. The modifications 
in both the CMB and matter power spectrum is originated thanks to the fact that, unlike in the cosmological 
constant case, the electromagnetic dark energy model produces fluctuations that might modify the evolution of the 
gravitational potential. In that figure, we can see that the small scales behavior is unaffected with respect to the 
standard ACDM case. The reason for this is that the electromagnetic perturbations decay very rapidly once they 
enter into the horizon as we have shown in the previous section so that only those electromagnetic modes whose scales 
have become sub-Hubble very recently (corresponding to the low multipoles part of the spectrum) can contribute in 
a non-negligible way to the metric perturbation evolution through Einstein equations. Moreover, since dark energy 
density is negligible during decoupling, the contribution of electromagnetic perturbations to the ordinary Sachs- Wolfe 
effect is also negligible. The main effect will appear in the late-time Integrated Sachs- Wolfe (ISW) effect as due 
to the evolution of the metric perturbation. The analytical estimate of such an effect is difficult to obtain since it 
requires to know the time evolution of the metric perturbation when the electromagnetic perturbations contribute in 
a non-negligible way to the Einstein equations, so that the approximated solutions obtained in the previous section 
are no longer valid. However, we can easily understand how the ISW effect will be modified by noticing that those 
modes that have become sub-Hubble very recently can still have an appreciable amplitude and, therefore, modify the 
late-time evolution of the gravitational potentials, which gives rise to a modification of the ISW, but only for low 
redshift so that the early ISW remains unaffected. Moreover, if the corresponding mode crosses the horizon with a 
too large amplitude, the modification in the gravitational potentials evolutions might be excessively large and, thus, 
conflict with observations. This is the reason why on large scales we obtain some distinctive signatures for large 
enough values of the primordial amplitude A. Notice that such signatures are more apparent in the matter power 
spectrum. Instead of giving the constraints in terms of the parameter A we shall give the results in terms of the more 
physical quantity 5 a defined as: 

-1/2 

Sa = ■ (34) 
PAo 



As we have already commented, in a more realistic quasi-de-Sitter inflationary epoch, the power spectrum of the temporal component 
becomes slightly red-tilted with a spectral index ua = 0(e) with e the slow-roll parameter }2gj. However, we shall neglect this small 
spectral index in order to set the initial conditions, being its effect expected to be small. 
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with Vk = \ Spk\ 2 the power spectrum of the electromagnetic energy density perturbations. Hence, the magnitude 
8a gives the amplitude of the energy density fluctuations of the electromagnetic field at a given scale k relative to 
the homogeneous contribution. Since 5pu and pa evolve in the same way on super-Hubble scales, the quantity 8a 
does not depend on time as long as the corresponding mode remains super-Hubble. Notice also that Spk contains 
two types of contributions, namely: one proportional to the metric perturbation and other proportional to the 
electromagnetic perturbation. Thus, in the case when the component proportional to the electromagnetic perturbation 
becomes dominant, the quantity 5a becomes scale-invariant on super-Hubble scales because Spk oc <L4ofc and 5Aok 
is proportional to fc~ 3 / 2 due to the flatness of its primordial power spectrum. However, if the metric perturbation 
contribution is dominant 5a will depend on the wave-number of the considered mode. In any case, the upper bound 
that we shall obtain for 5a will show how large the primordial electromagnetic perturbations arc allowed to be in order 
to be compatible with CMB measurements. In particular, we obtain that 5a must satisfy the constraint <5^~ 10~ 7 in 
order not to be in conflict with the CMB quadrupole. In fact, the overall effect on the CMB power spectrum would 
be that the higher value of 5a is, the more tilted the lower multipoles part of the spectrum becomes. In fact, since 
the 5pk modes decay rapidly once they reenter into the horizon (as we already commented above) the only important 
effect appears for the quadrupole and, as a consequence, the bound on 5a is actually a bound on such a quantity at 
the quadrupole scale. 

On the other hand, the upper bound obtained for 5a can be linked to a variation of the Hubble parameter in a 
quasi-de-Sitter inflationary epoch where the the Hubble parameter is not exactly constant but varies slightly. To show 
such a link, we have to remind that the background energy density of the electromagnetic field is given by all the 
modes whose scales are larger than the Hubble radius today, whereas the perturbations correspond to modes whose 
scales are smaller than the present Hubble radius, as commented before. In other words, the background is given by 
the modes that remain super-Hubble at the present epoch and the perturbations correspond to those modes which 
have already entered into the horizon. In [28| we showed that the amplitude of the electromagnetic field fluctuations 
(for the temporal component) at a scale k is given by Hk, with Hk the Hubble parameter at the time when the scale 
k exits from the horizon. Now, let us notice that 

Sa ~ Tpi — [wj (35) 

where Hi is the Hubble parameter at the beginning of inflation. Then, the constraint <5^4< 10 7 implies that the 
Hubble parameter must have reduced (at least) by a factor ~ 200 since the beginning of inflation until the time 
when the scale of the present Hubble radius (quadrupole scale) left the horizon, i.e., Hk ~ Hj/200. Since the Hubble 
parameter is proportional to the square of the inflation scale we also obtain that the scale of inflation must reduce 
in a factor ~ 15. Notice also that this fact requires a red-tilted spectrum for inflation, i.e., the Hubble parameter 
must decay throughout the inflationary epoch. Moreover, since 5Aok ~ -Hfc (as we have already said) we have that 
the fluctuations of the temporal component must satisfy: 

^<5x 10- 3 . (36) 
Ao 

Then, we conclude that the electromagnetic dark energy model is compatible with current CMB and LSS unless very 
large electromagnetic initial perturbations are generated during inflation. 



VI. DISCUSSION 



In the quantization of the electromagnetic field, either in the covariant or in the path-integral formalism, one has 
to introduce a gauge-fixing term in order to eliminate unphysical modes. The form of the gauge fixing term is unique 
in the sense that it is the only one that does not require the introduction of additional scales in the action and leads 
to linear equations of motion, although one may even consider nonlinear terms, in which case Faddeev-Popov ghosts 
appear also in flat space-time quantum electrodynamics. It has been proposed that the familiar, linear gauge-fixing 
term might have physical consequences in curved spacetime. One then promotes the gauge-fixing term into physical 
term which is included in the fundamental theory and not introduced during the quantization procedure as usually. In 
fiat space, the usual predictions of electromagnetism are retained. However, a gravitational field may excite a scalar 
degree of freedom of the electromagnetic field, which now is not a gauge mode but has physical consequences. A 
remarkable result of this is that in the homogeneous and isotropic background the temporal component of the vector 
A^ assumes an evolution that conspires with the cosmological expansion to contribute as an effective cosmological 
constant. This constant is of the order of magnitude required by the present cosmological data given that the field 
was generated during inflation occurring at the electroweak scale. 




FIG. 2: In this figure we show the CMB power spectrum (left panel) and the matter power spectrum (right panel) for both the 
electromagnetic dark energy model (dotted lines) and the standard ACDM model (solid lines). We have plotted several cases 
with increasing values of the initial amplitude for the electromagnetic perturbations and we see that the only modifications 
appear for large scales. In particular, the CMB quadrupole becomes excessive large, being incompatible with WMAP data 
(green dots) for 5a > 10~ 7 and the small k region of the matter power spectrum becomes very different from that of ACDM. 
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FIG. 3: In this figure we show the evolution of the gravitational potentials for small (k — 5.8 x 10~ s Mpc -1 ) and large 
(k — 1.1 x 10 -5 Mpc -1 ) scales in the left and right panels respectively. We also show the evolution in a ACDM model (solid 
blue lines) for comparison. As commented in the main text, the evolution for the small scales is exactly the same as in the 
case of a cosmological constant. In the large scales case we see how the evolution can be very different for large enough initial 
amplitudes of the electromagnetic perturbations. We plot the electromagnetic model in dotted (red) lines and we see that the 
larger the electromagnetic perturbations are, the more distinctive the evolution of the gravitational potential is. 



Such a term, being not covariantly constant but having internal dynamics, implies there will be fluctuations. We 
analysed these fluctuations both analytically and numerically. At superhorizon scales, the perturbations tend to freeze 
to constant value, while after entering horizon, they decay exhibiting oscillations. As the cosmological constant begins 
to contribute more significantly to the energy density, the fluctuations of the field begin to affect the gravitational 
potentials. Since this happens at recent times, the effects to the matter power spectrum is confined to large scales, 
and to the CMB at small multipoles. However, depending on the initial amplitude of the perturbations, which in turn 
depends on the details of inflation, these effects might be detectable. In terms of the parameter we find 5a ~ 10 -7 . 
This implies that the amplitude of fluctuations generated for the electromagnetic field is < 5 x 10~ 3 . Moreover, this 
fact predicts a reduction of the scale of inflation in a factor (at least) ~ 15 since the beginning of inflation until the 
scale of the quadrupole crossed the horizon. Since the mode V^A^ behaves as a free test scalar field, see Eq.©, we 
expect that it acquires the spectrum whose amplitude is given by the Hubble parameter just like that of the inflaton 
field (with just a slightly different slope due to negligible backreaction) . 

This shows that the model can be perfectly viable, but even more interestingly, for some specific inflationary 
scenarios there could be observable signatures from the model. Initial fluctuations of the order > 10 -3 can already 
be ruled out. In addition, if a coupling existed between the inflaton and the electromagnetic field, or if the temporal 
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component of the latter contributed non-negligibly to the inflationary expansion, one would expect their spectra to 
be correlated, which could yet lead to a new kind of signatures. 
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APPENDIX: PERTURBATION EQUATIONS IN THE SYNCHRONOUS GAUGE 



In this appendix we shall derive the equations of the perturbations and all the other expressions in the synchronous 
gauge, for which the perturbed line element is: 

ds 2 = a(r) 2 [dr 2 - + h^dxUx 1 ] (37) 

The scalar modes of the perturbation are usually expressed in terms of the scalar functions h and rj defined by means 
of HI: 



= I d A k 



fci ^(r)+f#-^i]<,M(» 



1/2 "'■ " ' ' \ j,i :->'"'' ' ' ' e%hX ( 38 ) 

To go from one gauge to another we perform a coordinate transformation given by: 

x" = x" + C 1 (39) 

As we are considering only scalar perturbations we can set ^ = (a, V/3). The relation between the scalar metric 
perturbations in both gauges is as follows [Hj]: 

4>k = ^K + 6r,^ + H(h' k + 6n' k )] (40) 

<t>k = m - + 6r?D (41) 



with a k = j3' k and 



Pk = -^(h k + 6r] k ) (42) 



The transformation for the vector field and its energy-momentum tensor are given by S^A^ = —CqA^ and 
S^T^l — —CcT^ u respectively. From these transformation laws we obtain for the vector field 

<5 c .4o = -K)' (43) 
5 C A % = -d.aAa (44) 

so that the perturbations SAq and SA in both gauges relate as: 

SA c Q onf = 5A s yn + (aA o y - 2HaA 
A conf _ A syn +aAQ ( 45 ) 

These relations show that the combination <5.4o — A + 2H.Ao is gauge-invariant (in the sense of cosmological pertur- 
bations). 

On the other hand, for the energy-momentum tensor we obtain: 

<5 C T° = -aT% (46) 
kT\ = -aT\ (47) 

= -(r%-±T^d ia (48) 

It is very interesting to note that, given that the background evolution of the vector field is the same as that of a 
cosmological constant we have that both the unperturbed energy density and pressure are constant and they satisfy 
p + p = so that the energy momentum tensor remains invariant after the gauge transformation. In fact, this is 
what one would expect from Stewart- Walker lemma [36| . Notice also that due to the gauge invariance of the energy 
momentum tensor the results showed in Fig. [T] are also valid for the synchronous gauge. 
The perturbed equations for the vector field in the synchronous gauge are: 

SA'^ k - mSA' Qk - 3k 2 dA 0k = -2k 2 (2SA' k - 5H6A k ) - - [(h' k ' - 4Hti k )A + h' k A' ] (49) 

SA' k ' - 4HSA' k + Uh 2 - 2H' - hA SA k = 2 (JsA' ok - HSA ok ) + ~h' k A (50) 
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and the perturbed energy-momentum components in the synchronous gauge are: 

ST° = ^{(25A' ok + A h k )A' () + 2k 2 {(-3SAo k + 35A k -6H6A k )An + A' 6A k }} 
ST) = ^ { {2SA' ok + A h' k ) A' + 2k 2 [3SAo k - 36A' k + 6HSA k )A Q + A' SA k ] } 8% 
5T% = ~^A [h' k A + 25A! ok + 2k 2 SA k ] (51) 
Condition (|25ll remains the same in this case: 

2Ank 2 

S(Pk + Pk) = Y~ ( SA ™ - SA 'k + 2WM fc ) (52) 

a 

which is just a consequence of the aforementioned fact that the combination 5Ao — 5 A' + 2HSA does not depend on 
the gauge choice. 



